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Recently, a series of noncentrosymmetric superconductors has been a subject of considerable 
interest since the discovery of superconductivity in CePtaSi. In noncentrosymmetric materials, the 
degeneracy of bands is lifted in the presence of spin-orbit coupling. This will bring about new effects 
in the Kondo effect since the band degeneracy plays an important role in the scattering of electrons 
by localized spins. We investigate the single-impurity Kondo problem in the presence of spin-orbit 
coupling. We examine the effect of spin-orbit coupling on the scattering of conduction electrons, by 
using the Green's function method, for the s-d Hamiltonian, with employing a decoupling procedure. 
As a result, we obtain a closed system of equations of Green's functions, from which we can calculate 
physical quantities. The Kondo temperature Tk is estimated from a singularity of Green's functions. 
We show that Tk is reduced as the spin-orbit coupling constant a is increased. When 2qA:f is 
comparable to or greater than ksTKia = 0), Tk shows an abrupt decrease as a result of the 
band splitting. This suggests a Kondo collapse accompanied with a sharp decrease of Tk- The 
log T-dependence of the resistivity will be concealed by the spin-orbit interaction. 
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I. INTRODUCTION 

The Kondo effect has attracted many researchers since 
the discovery of the solution of the resistance minimum 1, 
[3] . The effect arises from the interactions between a sin- 
gle magnetic atom and the many electrons in a metal. 
Metals, when magnetic atoms are added, and rare earth 
compounds exhibit many interesting phenomena that are 
related to the Kondo effect. The spin-flip scattering of 
a conduction-electron spin by a localized impurity spin 
gives rise to a term proportional to In T in the resistivity. 

Superconductors without inversion symmetry have at- 
tracted much attention since the discovery of supercon- 
ductivity in CePtaSifl]. A group of noncentrosymmet- 
ric rare-earth compounds has been reported to exhibit 
superconductivity: for example, Li2Pt3B[l,Q, CelrSis, 
CeCoGeg, CelrGeg , and LaNiCai- The absence of 
spatial inversion yields the splitting of bands due to a 
spin-orbit interactionfiol [T]| . 

The influence of the spin-orbit interaction was dis- 
cussed very recently in two-dimensional systems starting 
from the single-impurity Anderson model 12 14]. In the 
conventional Kondo problem, the conduction-electron 
states with spin up and down are degenerate. We expect 
that the band splitting has a large effect on the Kondo 
effect, and is closely related to a multi-channel Kondo 
problem. The purpose of this paper is to investigate this 
subject on the basis of the s-d Hamiltonian with the spin- 
orbit interaction of Rashba type in three dimensions at 
finite temperature. We calculate Green's functions and 
evaluate the Kondo temperature Tk from a singularity 
of them. We show that Tk is reduced as a result of the 
band splitting and shows a abrupt decrease when akp is 
comparable to ksTK- 

The paper is organized as follows. In Section II we 
show the Hamiltonian, and in Section HI we derive equa- 
tions for Green's functions. We obtain an approximate 
solution in Section IV. The Kondo temperature and cor- 



rection to resistivity are discussed in subsequent Sections 
V and VI. In Section VII we examine the strong limit of 
the spin-orbit interaction where the details of calcuations 
are shown in Appexdix. 



II. MODEL HAMILTONIAN 



The Hamiltonian is H = Hq + Hsd = Hk + Hgo + Hsd 
where 



Hk 



Xl?k(4tCkT + 4iCki), (1) 

k 

^[a(?fca; -I- ky)c\^^c\ci + a{-ikx + A:y)cJ^j,Ckt], 

k 

(2) 



2 N 



't 



kk' 



-S-cl^Ck'i]. 



(3) 



is defined by ^k = Ck ^ M where ek is the dispersion 
relation of the conduction electrons and fi is the chemi- 
cal potential. Cko- and are annihilation and creation 
operators, respectively. We set Hq = Hk + Hgo- S+, 
S- and Sz denote the operators of the localized spin. 
We consider the spin-orbit interaction of Rashba type in 
Hso- OL indicates the coupling constant of the spin-orbit 
interaction. The term Had indicates the s-d interaction 
between the conduction electrons and the localized spin, 
with the coupling constant J. J is negative for the anti- 
ferromagnetic interaction. 
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III. GREEN'S FUNCTIONS 

First, we define Green's functions of the conduction 
electrons 

Gkk'.(r) = -(r.Ck.(T)4^(0)), (4) 
i^kk'(r) - -(T.ck;(r)4,^(0)), (5) 

where T^- is the time ordering operator. We note that the 
spin operators satisfy the following relations: 

S±S:, = T^Sz, SzS± = ±^S±, (6) 
S+S- = ^ + S,-Sl (7) 
S-S+ = l-S,-Sl (8) 

We also define Green's functions which include the lo- 
calized spins as well as the conduction electron opera- 
tors. They are for example, following the notation of 
Zubarev|15|, 

((5.Ckt;4't))- = -(r.5,Ckt(r)4,^(0)), (9) 

((5_cki;4^)). = -(r.5_Ck^(T)4,^(0)), (10) 

{{S.c^i;cl^})r = -(T,5,cki(r)4,^(0)), (11) 

((5_Ckt;4,^)). = -(rr5_Ckt(r)4,^(0)). (12) 

The Fourier transforms are defined as usual: 

Gkk'.(r) = i^e-^""^Gkk'.(«c<.„), (13) 

n 

i^kk'(T) = ^5]e-^""-i^kk'(*^„), (14) 

((5.ckt;4,t))r = ^E^"""'«^-^kt;cLt)).-Ji5) 

n 

(16) 



From the commutation relations 



Similarly, the equation of motion for -Fkk' is 
d 



dr 



Fkk'ir) = -Ck^kk' - a{-ikj: + fcy)Gkk't('^) 



(20) 

We define 

rkk'(r) - ^^e-'""rkk'(«^„) 

n 

= ((5.ckt;4't))- + ((^-Cki;4,^))., (21) 

^^e— "-$qk(*a.„) 



'i'qk(T) 



= {{SzCqi - S'+Cqt;4.^)), 



(22) 



then we obtain 



(it^ri - Ck)Gkk't(*'^n) = 5kk' + "(^^a; + fcy)^kk'(»W„) 

- -^5]rqk'(»c.„), (23) 



qk' 



- ^k)^kk'(*'*^n) = a(-«fcx + fcy)Gkk't(«^n) 

+ ^E'^qk'(«^„). (24) 

9 

To obtain the solution to the equations above, we need 
Green's functions in eqs.(9)-(12). The equations of mo- 
tion for ((S'jCkt; C|{^,^)) and ((S'-Ckj.; cj^.,^)) are 

{iu - Ck)((S'2Ckt;4't))'" 

- j^Yl [(('S'fcqt;4't))'''^ + ^(('5-Cqi;4't))* 

q 



[Ho, Ckf] = -CkCkf - a{ikx + fciy)ck;, (17) 
[Hsd,cwt] = -^Xl^^-^^^k't - 5'-Ck'i), (18) 



the equation of motion for Gkk't(''') reads 
d 

^Gkk't('r) = -<5(T)(5kk' - CkGkk't(T) 
-a{ikx + fcy)-Fkk'(T) 



{{S-Cktcl^^Cq>l;cl^,^))^ 



{iuj - Ck)((5'-Ck;;4't))*'^ 
a{~ik.^ + ky){{S-C]^^;cl,^))i 



(25) 



qq' 



[(('S'-Cki4t^q't;4't))» 



qq' 



((S'_Ck^4^Cq-;;4,^))i„ - 2((S'2Ck;4j^Cq't;4,^)) 



•qi^q T' ^k't'/*" 



2iV 



(19) 



(26) 



q' 
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We use the commutation relation S+S- = 3/4 + S'z — S'f 
to obtain 



{ioj - Ck)rkk'(«w) 

= Sww'{Sz) + a{ikx + fcy)((5'2Ck^;ci.,^)) 

J r3 



~ [(('^+Ckt4iCq'T;cLt))i 



qq' 

— ((S'-Ck-fCq^Cq'J,; cjj.,^))itj — ( (5'^ Ck4,Cq^Cq'-f-; cj^.,^)) itj 
+ ((S'-Ck4c|j^Cq'^;4,^))j^ 

+ 2((5,Cki4cq.t;4't))-]- (27) 

Here we assume that {Sz) — 0. Now we need 
((S^Ckj,; cj^/j-)) and ((5-Ckt; cjj.,|)) to obtain a solu- 
tion for Fkk'- The equations for ((S'^Ckj,; cj^/-|~)) and 
((5-Ckt; cj^/|)) read 

(iw - Ck)((S'zCki;4,^)),,^ 
= + fc;,)((5^ck_i.;cj^,|))ic^ 

qq' 

- {{S-c\^Cc^'iCl,Y,c]^,^))^u^ , (28) 

[iuj - ^k)((5'-Ckt;4't))»" 
= -5kk'(<S'-) +a(«fc^ + A;y)((5'_ck4,;4,|))iaj 

+ ^II(('^-Vt;Ck't))^'^ 
q' 

qq' 

- ((5'-Ckt4^Cq't; c['t))«" + 2((5'24j_cq-^ckt; c[,^)), 

(29) 



IV. APPROXIMATE SOLUTION 

We assume that the spin-orbit coupling a is small and 
we keep terms up to the order of a. We adopt the ap- 



proximation that 



«w - fk 



((^2Ckt;4,^))i 



(30) 

((5_ckt;4,,)).. = ^^^^^-i^((5-ck;;4,,))...(3l) 

«w — 4k 



This means that we have neglected the terms of the order 
of J a in the right-hand side. Then we obtain 



iuj - ^k - 



'rkk'(ia;) 



IUJ - 4k 
[^((cqt;4't))"^ +rqk'(ja;) 



qq' 

- ((S'_Ckt4^Cq'4,; 4't>>''" - ((5'-Cki4tCq't; 4't))'" 
+ (('S-Ckiclj^Cq'^; 4't))«" + 2((52Cki4^Cq.^; 4't»»'^ 

(32) 

We use the same approximation in the right-hand side 
of eq.(23), that is, {iuj — Ck)-Fkk'(*w) = a{—ikx + 
ky)G]i]i'^(iuj), and we have 

{iuj - ^k)Gkk't(«^) = Sww " ^ X! rpk'(jw) 

p 

+ ^ " , "^ Gkk't(»^) (33) 
lOJ - 4k 

Here, we employ the decoupling approximation proce- 
dure for Green's functions (la. \IT\ . Many-body Green's 
functions are approximated as follows. 

((5'-Ckt4tCq'i;4't))''" ~ (Ckt4t)(('5'-Cq.^;4'^))jc^ 
+ ('S'-4^Cq'4,)((Ckt;4't))»"' 

(34) 

((S'+ckt4^Cq't;4't))''" ~ (S'+4^Cq/^)((ckt;4't))«'^ 

- ('5+4iCk't)((Cq't;4't))»'^' 

(35) 

((5'^Cki4^Cq't;4/t))''" ~ (cki4^)((5'^Cq't;4't))''" 

- (5^4^Cki)((Cq't;4't))«"' 

(36) 

((S'_Cki4j^Cq'4,;4't))''" ~ (Cqi4'i)(('S'-Cki;4't))'" 
+ (CkiCqi)(('5'-Cq'i;4't))'"- 

(37) 
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We define 

q q 

ruk = 3^(5_4cH)=2^((5,4ckt) + (5-4cH)). 
q q 

(39) 

We used the relation obtained from the rotational sym- 
metry in the spin space, 

(40) 

Then, after the analytic continuation — >■ w + iS, we 
have 

^ ^II^qk'(w)+ Uk- -^^rqk'(w) = 0, 

q ^ ' q 

(41) 



we obtain 



(38) rkk'(a;) = —Gi{uj)Gi, {u) ( mk - ^ ) (1 + JG{cj)) 



(49) 



«k - 2 ) -^r(a;) 



1 + JG(w) + ( J/2)2r(w)F(w) ' 
Gkk'(a;) = 5kk'G0(u;)-^r(a;)G0(a;)G0,(a;) 



1 + JG(w) + ( J/2)2r(cj)i^(w) 
5kk'G^(u;) + ^G^(c.)G^.(c^)i(a;), (50) 



where we defined 
t{uj) = 



r(^) 



4 1 + JG{lu) + ( J/2)2r(w)F(w) " 



(51) 



rrik is given by 



2E((^-4tCqt)(^-4tCq^)) = 2Erqk(T = "0) 



= ^Ee""'l'qk(ia;„), 



(52) 



(42) 

where we set fc^ = k'^+ky. Then, we obtain from eqs.(41) 
and (42) that 

rkk'(a;) = GO(a;)(mk-J) ^^k'M 

- G° (c.) [Lu-l)j+ (mk - ^) (^) ' ^('^) 



where 5 is an infinitesimal constant. Because mk is real, 
we obtain 

mk = -4^Ee''^"*G0(za;„)i(*'^n). (53) 

Similarly we have 

rik = EG^(ia;„)(l + JF(iw„)t(ia;„)). (54) 



(43) 



where 



V. KONDO TEMPERATURE 

A singularity of t{uj) determines the characteristic 
temperature of the system. We investigate the high- 
temperature region where rrik = 0. Then, 



2 \w — ^k + oik± w — ^k — cek± 



nk 



G°(a;), 



(44) 
(45) 

(46) 

(47) 



r(a;) = --F(w). 



(55) 



We obtain for k = k' 
Gkk(w)-' = Gi{io) 



_i 3J2 F{o 



16N 1 + JG{ij) 



(56) 

The Kondo temperature is determined by the van- 
ishing of the denominator in this equation: 



Because of 



+ 



w — ^k + o:k± w — ^k — ak± 



^ErqkM-^r(-)^^MTT7GHT(J 



{J/2rT{c.)Fiu;)' - 4V 
(48) 



1 ftanhf^!!^ 



2rj 



if 



tanh 



^k + ak_ 



2T, 



K 



= 0, 
(57) 
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where 



(58) 



This yields the temperature Tk as 

2e^D 



ksTK 



■ exp 



We have used rik = (/(^fc - aki_) + /(^fe + ak±))/2 by 
neglecting the term of the order of J^. By using the 
expansion, 



1 _ 7C(3) / 2akF 
Pf\J\ ~^^\kBTK 



31C(5) (2akF 
ksTx 



127r4 



(65) 



tanh 



z — m{2n + 1) 



(59) 



where we introduced the Boltzmann constant ks and 
the density of states pp = rnkp /{2Tr'^). This is a self- 
consistency equation for Tk, and yields 



the equation for Tk is 



1 



8(27r) 



K 



/ dk±k_i 



+ 
+ 



M Jo 

ZTrT/f sign(2n + 1) 
to + 2ak± — iTr(2n + 1)Tk 
i'7rTKsign{2n + 1) 



«7rTxsign(2n + 1) 



exp 



- 0.21314 



fa,. 

0.0550 — 
V X 
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uj-m{2n+l)TK - 0.01655 ( — ) + 0.005396 (—) - 0.001822 ( — ) 



LJ - 2akj_ - m{2n + 1)Tk 
where ii" is a cutoff and we use an approximate expression 



(60) 



f 

J-K 



K LO - kj_/{2m) + ak_L - k^/{2m) + iJ. 



k±/{2m) + hll{2m) + akj_- ii- m{2n + 1)T 



K 

2m 1 

i7rsign(2n + 1)Tk ; 7 zt^- 

H ^ ^ ^ ' uj + 2ak±-iTr{2n+l)TK 

(61) 

Wc set an cutoff no = D/ {2'kTk) in the summation with 
respect to n. By using the formula for the digamma 
function, 



1 



71=0 ^ 



•0 ( ^ + a; + no 



V ( ^ +a; 



(62) 



we obtain 
1 = |J| 



f2eiD\ (\ 



327r2 

1/1 a; + 2aki_ 
2^ yi i2TxTK 

\^ (\ + ^2:^^ -\^(\- ""rjf^ ) 



i2-KTK 



i2iTTK 



I I (I w - 2afcj 
2J V 2 i2iTTK 



(63) 



We set fi = kp/{2m) and K = 2kp, and expand the 
digamma function in terms of ak±/ (2TrTK)- For w = 0, 
we have 



\J 



mkp 



loa 



2eiD 
■kTk 



7C(3) (2ak 



27r2 V Tk 



+ 



31C(5) (2akF 



127r4 V Tk 



(64) 



+ •• 

= 9{x), 
with variables 



(66) 



Tk/TI 



lO 



— 2akp/kBTK, 



where 



ksTK = exp 



Pf\J\ 



(67) 



(68) 



We expanded g{x) in powers of ar/x up to the tenth 
order. The function g{x) is shown in Fig.l, where higher- 
order terms are small and negligible except near x ~ 
0. The equation x = g{x) has no finite solution when 
ar > 1.045. This indicates that Tk vanishes when the 
spin-orbit coupling akp is greater than 1.045A'bT]|., and 
indicates a Kondo collapse with a sharp decrease of Tk- 
This may overestimate the reduction of Tk- When a is 
very large, if we use the asymptotic relation ^(1/2 + 2) ~ 
log(z), we obtain 



1^Pf\J\ 



logf^)-ilogf^Ui 
^\nkBTj 2 ^yTTkBTj 4 



(69) 



This yields 



ksTi 



K 



{2e'^Df 
2akF TT 



■ exp - 



Pf\J\ 



(70) 

for ar ^ 1- We show Tk as a function of ar in Fig. 2. 

We expect that, in the strong coupling limit ar ^ 1, 
Tk should approach that of single-band model: 



ksT^ 



2eiD 



■ exp 



Pf\J\ 



(71) 



We will show this in the section 7. This agrees with 
eq.(70) for akp ~ D. This is very small compared to 
the original Tk because T^/T% ~ ksT^/D. Therefore 
Tk decreases as is increased and shows up a sharp 
decrease at ~ 1- 



6 



1.5 



1 - 



g 



0.5 







n — '~i 




r ^ 




/ A 1 




: 1 J \ J 





0.5 



1.5 



X 



FIG. 1: g{x) = exp(-(7C(3)/27r2)(a^/a;)2 + • • • ) up to the 
tenth order of Or/x as a function of x for Or = 0.5, 1 and 1.5. 
The straight line is a hnear function x. 




FIG. 2: Tk as a function of Qr. x = has no solution for 
Or > 1.045 (Kendo collapse) as indicated by dotted line. The 
dashed line is an expected line. 



VI. CORRECTION TO RESISTIVITY 



The imaginary part of G'kk('^) ^ gives the scattering 
rate of conduction electrons due to the localized spin. 



The inverse of the life time Tk{oj) is 
1 



ni7VImGkk(w)"^ 



16 



1 + JK{uj) 



{l + JK{oj)y + {JL{uj))^' 

(72) 



where rii is the concentration of magnetic impurities. We 
have defined K{oj) = ^eG{uj+i5), L{lo) = —lm.G{uj+i6), 
and p°^{oci) = — (l/7r)ImF(a; + iS). Because we obtain 



, f 2e<D\ 7C(3) f2akF 



47r2 \kBT 
the formula of the conductivity yields 

O- = y J Tfe(Cfe)'Uj.— p(Cfe)d^fe 

16 



(73) 



2e^ 
3 

2e2 



t^Fp(O) 



37mjJV"(0) 
3 ^^^37rn,|J|/5"(0) 



(1 - \J\K{0)) 



T 



+ 



7C(3) /^2afcF 
47r2 



(74) 



We have the term (a/T)^ that comes from the spin-orbit 
interaction, and this term will conceal the logarithmic 

dependence of the resistivity. Then, the electrical resis- 
tivity R in the high temperature region T ^ is 



R 



Ro 



l-\J\p^log(^] 



(75) 



where Rq is a constant. If the term (a/kBT)'^ is larger 
than the logarithmic term, the resistivity even shows 
R ^ T^. Hence, the spin-orbit coupling may change the 
temperature dependence of the resistivity drastically. 



VII. STRONG SPIN-ORBIT COUPLING CASE 

In this section let us consider the case with strong 
spin-orbit coupling. For this purpose, we diagonalize the 
Hamiltonian Hq: 



Ho - E(4tcU) ( ai-tt+ ky) a J \ Cki 



a{ikx + ky) \ f Ckt 
^ [(a - afc±)4ak + (a + afc±)/3^/3k] , (76) 



where k±_ = \ k'^ + k^, and ak and /3k are defined by 



Ok = UkCkt + "kCki, 
/3k = -^'kCkt + UkCkj,. 



(77) 
(78) 
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The coefficients and Wk are 



VIII. DISCUSSION 



(79) 



satisfying bkP = 1- We consider the case where the 
band split is so large that we can neglect the upper band. 
This means that we keep terms that contain a-operators 
only. In this approximation the s-d interaction term is 



J 1 



kk' 



(80) 



This is the model of one-channel conduction-electron 
band that interacts with the localized spin. 
Let us consider the following Green's function: 



G£k'M--(T.ak(r)4,(0)). 



(81) 



By using the same method in previous sections, Ggj,, is 
shown to be 



<5kk' 



J 



i'ki'k' 



2: - CkQ 2iV (z - ^ka)(z - Ck'a) 



for arbitrary complex number z where we defined 
3 J F„(z) 



t{z) 



16l + ^G„(z)-^(^)V„(2)2' 



^ ly"k.-l/2 



(82) 
(83) 

(84) 
(85) 



We derive this formula in Appendix. The Kondo tem- 
perature is determined from a singularity of t{z) in 
the same way as previous sections. We obtain 



We investigated the Kondo effect in the presence of 
spin-orbit coupling. The influence of band splitting was 
examined by using the Green's function method where we 
adopted the decoupling scheme to obtain an approximate 
solution. The Kondo temperature is reduced by the spin- 
orbit interaction, and shows a sudden decrease when akp 
is of the order of ksT^. We call this the Kondo collapse 
due to the spin-orbit coupling. The Kondo effect is sup- 
pressed and the log T-dependence of the resistivity will 
be weakened and concealed. The reduction of Tk as a re- 
sult of the spin-orbit coupling is consistent with the result 
for the single-impurity Anderson model using the numer- 
ical renormalization group technique [13|. In their work 
the Kondo temperature is a decreasing function of the 
Rashba energy En oc when the level of the localized 
electrons is lowered, that is, in the Kondo region, while 
it is a increasing function when the localized level is not 
deep. The variation of the Kondo temperature is approx- 
imately linear as a function of En, namely, quadratic in 
a. This is consistent with our result which shows a small 
variation of the Kondo temperature with the quadratic 
correction when a is small. In a recent work[l^. the 
Kondo temperature is increased in the presence of the 
Dzyaloshinski-Moriya interaction. The Dzyaloshinsky- 
Moriya interaction, however, vanishes in the Kondo re- 
gion with particle- hole symmetry = —U/2. Hence the 
result in ref. 121 seems consistent with the result for the 



s-d model. 

As a limit of strong spin-orbit interaction, we can inves- 
tigate a crossover to a one-channel Kondo problem. The 
Kondo problem with the spin-orbit coupling is closely 
related to a multi-channel Kondo problem. The Kondo 
temperature is reduced considerably because the degen- 
eracy of the conducting electrons becomes half of the 
conventional Kondo system in this limit. The specific 
heat also exhibits a log T-term in the present model with 
one-channel conduction band, and this term appears in 
the fifth-order of pj. This agrees with the original Kondo 
problem. 

The author expresses his sincere thanks to K. Yamaji, 
I. Hase and J. Kondo for helpful discussion. 



2eT D 
ksT^ = exp 



\J\PF 



(86) 



The characteristic energy is reduced significantly 
compared to the conventional Kondo temperature by fac- 
tor 2 in the exponential function. This factor appears be- 
cause the number of channel of the conduction electrons 
in this case is just half of the normal Kondo system. The 
resistivity is also calculated as 



IX. APPENDIX 

In this appendix we derive the equation of motion for 
Green's functions for the single-band s-d model and dis- 
cuss its physical properties. 

A. Green's functions 



R — -Ro 



with a factor 1/2. 



2 ^\7TkBT 



(87) 



Hq was diagonalized by ak and /3k : 

ttk = UkCkt + VkCkl, 

/3k = -I'kCkt + "kCkj,- 



(89) 
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The coefficients Uk and are 



1 'ikx ~\~ ky 

■Wk=;r, -^k = — 



V2k± ' 



(90) 



satisfying Uk + kk| = 1- The s-d interaction part be- 
comes 



kk' 



- {ukUk' - f^fei^fc')/5k/3k' 

- {ukVk' + VkUk')al0k' - {ukVy + UfeUfe')/3j^Q!k'} 
+ S+{vkUk'al.ai,' -UkVk'l^lPk' -VkVk'all3k' 

+ UkUk'l3l^au') 

+ S-{ukvl,alai,' - v^Uk'^lPw + WfeUfe'aj^/3k' 

- vlvk'plau')]. (91) 

The single-band s-d model contains only the following s-d 
interaction, 

Hsd = -^j^Yl^S;,{{ukUk' -Vkvl,)alak' 

kk' 

+ S+VkUk'olak' + S-Ukvl,alak']- (92) 

We consider the following Green's functions: 

G£k'(T) = -(r.ak(T)4,(0)), (93) 
((5,ak;4'))r = -(T,(5,ak)(T)4,(0)), (94) 
{{S+au;al))r = -(T,(5+ak)(r)4,(0)). (95) 

The Fourier transforms are defined similarly: 

G£k'(T) = ^^e--G£k'(z.;„), (96) 

((5.ak;4'))r = ^5^e'-"-((5.ak;4'))a., (97) 

((5+ak;4'))r = -^5]e-"-((5+ak;4,))c.. (98) 



The equations of motion for these Green's functions are 
derived as follows. 



^'^n<^kk'(«'^n) = ^kk' +^kaGkk'(iw„) 

~ ^Zl[("fe"9-'yfcO(('S'^aq;4'))i 



+ VkUq{{S+aq;a'l,))i,^^ 
+ UkV*{{S-aq;al,))iu,J, 



«w„((S'^ak; 4'))»w„ = ^kc<(('S'zak;4'))iwn 



+ VkU*g{{S+{nka- ^)aq,al,))iu^ 



(100) 



ia;„((5+ak;4'))»w„ = ^kc«((5'+ak; 4'))»w„ 
q 

+ 2ufeU*((SznkaQ:q;4'))iu-n " ^^fc^q (('S'+5'-aq; 4'))jc^„] > 

(101) 



ia;„((S_ak;4'))'w„ = ^ka(('S'+Q!k; 4'))iw„ 

+ ^X![("*^"9~^'=^9)^^'^-("ka - ^)Q;q;4'»»'^,. 

q 

- '^VkUg{{Szrn^aaq,Oit'))iujn - VkUq{{S-S+aci;al,))i^J, 

(102) 

where nka = ak^k- We have unknown functions 
((5'ankaQ!q;4')) a = z, + and -. 



B. Approximate Solution 

To obtain a consistent solution, we adopt the following 
approximation: 

((5'ankaQ!q;4')) = ("•ka) ( (-SaQIq; 4' )) • (1^3) 

Using this approximation and the relation S+S- = 3/4+ 
Sz — Si, we obtain 



(«W„ - ^ka)G£k'(«^«) = <5i 



kk' 



(99) 



+ MfcW*/((S'-aq';4'))i"„ 

+ (wfeUg' - ffcU*')((5'zaq';4'))»Wr.] 

+ ^(2^) E— ^Bw+^KOG^'k'i^'^.), 

(104) 
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where riq = (nqa). Here we define 
q 



Then Ggj^, and Fkk' read 

Wk'(2) = — + 



2;-^ka 2Ar(2;-^ka)(^-^k'a) 



t(z), 



(105) rkk'(^) 



This quantity reads after substituting the equations for 
Fkk' - ia;„ - xika ^ [ 

q q' 

- UkV*, (nq- i ) ((5'_aq';a[,))ia;„ 



(112) 
(113) 



for arbitrary complex number z. The Kondo temperature 
is determined from a singularity of t{z) in the same 
way as previous sections. We obtain 



Tia 
K — 



2e^D 



■ exp 



\J\pfJ ■ 



(114) 



The characteristic energy is reduced significantly 
compared to the conventional Kondo temperature by fac- 
tor 2 in the exponential function: 



- (ufev - VkV*,) { nq - ^ ) {{S^ac^r,al,)) 



[ukUq' + t;feW*/)G'q/k'(ia;„) 



(115) 



J 3 n-k — 1/2 
2iV8^ zw„-5q„ 



O.N 



This factor appears because the number of channel of the 
conduction electrons in this case is just half of the normal 



2-^ 8 ^ iujn — Koudo system. The resistivity is also calculated as 



X Y^{ukUg' +VkV*,)G'^,^,{i^0n)■ 
1' 

Then we obtain 



(106) 



R = Rq 



1 + 



Pf\J\, f2e^D\ 



log 



[nkBTJ 



+ ■■■ 



(116) 



with a factor 1/2. 



3 / J 



ii^n-^ka S\2NJ iuJn-^Ma 

X V I I 

^ia;„-^q,.l + _^2^|^ 

X Y^{UkUg + VkVl)Gl^,{iUn). (107) 

q 

We have set Wk = l/v^- Because z;k satisfies = —v-k 
and |ukp = 1/2, we have 



iOJn — Ck'. 



3 ( J\ 1 



8 V 2 



C. Entropy and Specific Heat 

The energy expectation value E = (H) is given by 

E = '^S.kialau} ~ ^^^ilS^iukUk' - Vkvl,) 

k kk' 

+ S+VkUk' + S^Ukvl,}ala]^') 

= ^Yl €kaGL(«'^n) -^^Yl rkk(iw„) 



2 pN 



iuJnt{iLOn) 



1 



l + (J/2)Ga(ia;„) 



where we set 



^ 2 - Cka ' 



We define 



i(^) = 



3J 



J^a(2) 



16 



l + ^G«(z)-^(^)^F„(z)2 



(108) 

(109) 
(110) 

(111) 



_ 1 \ ^ ^ka J 1 \ ^ tu.'n'-V 

(117) 

The expectation value of the interaction Hamiltonian is 
denoted as V is given by 

V = -■^Y.iiSziukUk' - VkVk') + S+VkUk' 



kk' 



+ S-Ukvl,}alau') 
J 1 



t{iuJri) 



2;3ArErkk(ia;„) = ^^^,,„_^,^' 

k kai„ 



(118) 
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This is written as 
J 



V = |p(0)Re / dujf{uj)t{uj - i5), 



where we adopted the approximation 
Fa{Lo ± i5) = =F7rp(0)«. 



(119) 



(120) 



p{uj) is the density of states of conduction electrons. 

We need t{z) to estimate V . Ga{z), which appears in 
the denominator of t(z), contains a singularity. Ga{z) is 
written as 



G„(z) - i?„(z) + -i-^^ 



1 



t{iuJn) 



2N PN ^ Z- ^ka {iuJn - Cka)^ 



(121) 



To estimate V", we use the expansion formula for the 
Fermi distribution function f{uj): 



du f {uj)h{Lj) 



-D 



dLoh{u) + —{kBTyh'iO), 
l-D 6 

(129) 

for a differentiable function h{Lu). Using this, we obtain 



V 



_37r 
~ 16 

37r TT^ 

16 Y 



p(0)JIm / du 



1 



{kBTyp{0)Jlni 



^ log{T^/T)-gi(3u;) 
d 1 



du; iog(r^/r) - ff(/3c^) 



=0 

(130) 



We are interested in logarithmic terms log(£'/fcBT), 
\og{D/kBTf and so on in the region | log(T^/T)| > 1. 
The second term is written as 



where 



1 



Fa{iuJn) - Fa{z) 1 



F,(z). (122) 



Ra{z) is evaluated as[l 
Ra{uJ-i5) ~ p(0) 



^ ^ 2 2tt 1 ^ I 2 27r 



(123) 

where ip is the digamma function and D is the cutoff 
energy. We use the following relation, 



1 



p(0)J 



, , D \ , /I .I3lj 
^°s(^^j-V'(^2+^2^ 



p(0)J 



(124) 



= -^fe^^p(0)'^ (iog(4/r))^ - 



(131) 



This is expanded as in terms of p{0)J: 



v. 



p{0)J 



( D \ 



TT 

32 V 2 y "° ^A^j 

37r^ /p(0)jy 
64 I 2 



(132) 



In the first term of V , the logarithmic corrections never 
emerge from the region where /3a; is large because we 
have g{j3u)) ~ log(/3u;) for large w and the T-dependence 
is canceled with log(Tj^/T). When f3u) is small, g{f3uj) is 
expressed in a power series of A dominant contri- 
bution is of the order of (log(T^/r))^^. The integral is 
restricted on the interval (— fc^T, 0) and the first term Vi 
is estimated as 



where 



(125) 



Then the interaction energy is 



V 



— p(0)JIm / duf{uj)- — -— 

16 J_D log(r^/r) - g{l3uj) 



(126) 

where we neglected the term of the order of (p(0) J)^ in 
the denominator of t{z). V has a logarithmic temper- 
ature dependence. Because of the relation between the 
free energy and V, 



V 



dF 

a/' 



(127) 



the additional entropy AS{T) is 



AS{T) - - — {F^Fo) ^- I — — V{J',T). (128) 



dJ' d 



J' dT 



1 



ksT 



J T , I 1 -^^ 

dujliaip - +1 — 
Z zvr 



-|^-^^(«)^(iog(Tij/r))^' 



0.00738 + 0.000026- 



1 



0.0052 



As a result, is given as 
A 



^ = -2^-^^(«)^(iog(T,?/r))- 



(133) 



(134) 



for a constant A > 0. 

From the relation — Dexp{2/{p{0)J)), we have 



dp{0)J 
7(0)7 



1 



logiT-/D) 



dlogT^. 



(135) 



Using this formula, the entropy obtained from the inter- 
action energy V is 



AS- 



dAF 



(136) 
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where 
AF = 



This suggests thatfl^ 



2 



^ 1 ( pmj 

{\og{D/T)f\ 2 



1 



log{D/T) 



log 



p{0)J 



log 



±K_ 

T 



iog(Tj^/r) 

(137) 



is the free energy. Because of the relation 



37r3 
' 4 
37r3 



(1 
1 



pJlog{D/kBTy 



4 (log(T;,/T))4 



(143) 



log 



±K_ 

T 



p(0)J 



log 



2e'f D 
n kfiT 



(138) 



up to the fifth order of p(0) J, AS* is given as 



AS = kuA 



1 / p(0)J 
3 V 2 

p(0)J 



2 

p(0)J 



1 f p{0)J 



1 

2 V 



log 
log 



D 
k^ 

D 
k^ 



3 |^p(o)jy 



log 



The logarithmic term first appears in the fourth order of 
p(0)J. Then the correction to the specific heat AC = 
TdAS/dT is 



as an expansion in terms of 1/ log(TR-/T). In the present 
model, the coefficients are reduced, where 4 is reduced to 
12/5 in front of pJ\og{D /ksT) in the specific heat com- 
pared to the usual s-d model, and the divergence near 
|he Kondo temperature is moderated. Because the for- 
mation of a local singlet by the the conduction electrons 
ksT J is weakened in a one-channel case, the entropy decreases 
more slowly as the temperature is decreased. 
(139) 

In the region | log(D/A:BT)| > 1 and |log(r^/r)| > 
1, AS" is obtained as a double-power series of 
ll\og{D/kBT) and l/log(r^/T): 



AC 

ks 



A 
2" 



p(0)J 



1 - 



12 / p(0)J 
5 V 2 



log 



D 
k^ 



(140) 

Hence the specific heat exhibits a logarithmic behavior at 
low temperatures. A logT-term appears in the fifth order 
of p{0)J; this agrees with the original Kondo problem^]. 
In the original Kondo problem, the entropy and the spe- 
cific heat were evaluated as[2l. [Toj 



AS^ksA 



1 



.{logiD /kBT)r 
Then we obtain 



log(Ti^/T) 



AC - fcR^ 



1 



1 



ASsd 
ACsd 



kB — {pJf 



3^^ 



{pjy 



3pJ log 
- 4pJlog 



D 
k^ 
D 



(141) 
(142) 



ilog{D/kBT))^ {\og{T^/T)) 



(144) 



(145) 
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